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Abstract 

Gravity is the most familiar force at our natural length scale. However, it 
is still exotic from the view point of particle physics. The first experimental 
study of quantum effects under gravity was performed using a cold neu¬ 
tron beam in 1975. Following this, an investigation of gravitationally bound 
quantum states using ultracold neutrons was started in 2002. This quantum 
bound system is now well understood, and one can use it as a tunable tool 
to probe gravity. In this paper, we review a recent measurement of position- 
space wave functions of such gravitationally bound states, and discuss issues 
related to this analysis, such as neutron loss models in a thin neutron guide, 
the formulation of phase space quantum mechanics, and UCN position sensi¬ 
tive detectors. The quantum modulation of neutron bound states measured 
in this experiment shows good agreement with the prediction from quantum 
mechanics. 

Introduction 

Phenomena due to the gravitational held have been well understood at the 
macroscopic scales. However, there are only a few cases of experiments at 
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microscopic scales, due to gravity’s extreme weakness when compared to the 
other forces, such as electromagnetic and nuclear forces. 

Even though quantum mechanics was established in the early 1900s, the 
Erst experiment to investigate a quantum effect under gravity was reported in 
1975 by the group of R. Colclla jT], in where a neutron interference pattern 
induced by a gravitational potential was observed. The major systematic 
uncertainty was attributed to a bending effect of an interferometer when one 
rotates the system to manipulate a relative quantum phase difference be¬ 
tween two neutron paths. A subsequent series of experiments were carried 
out with careful studies of systematic effects by measuring the bending effect 
using X-rays and considering the Sagnac Effect [2] due to the rotation of 
the Earth 011] • The latest attempt was performed using normal-symmetric 
and skew-symmetric interferometers and a two-wavelength difference mea¬ 
surement method. The agreement with theoretical calculation is at the 1% 
level [5j. 

Another approach to observe quantum effects under gravity was carried 
out by the group of V. V. Nesvizhevsky 00- In this experiment, one mea¬ 
sured a transmission of ultracold neutrons (UCNs) through a vertically thin 
neutron guide in the terrestrial gravitational held as a function of the slit 
thickness, and showed the evidence of quantum bound states in the gravita¬ 
tional potential by analyzing the minimum height of a guide through which 
UCN can be transmitted. In a later detailed analysis, the characteristic sizes 
of position-space wave functions for the first and second quantum states were 
evaluated with an uncertainty of around 10% [8j. The major systematic un¬ 
certainty in this evaluation came from the difficulty in modeling the rough 
surface of a scatterer used on the ceiling of the guide, which causes some 
deformation of the wave functions. To minimize the effects of the deforma¬ 
tion, an “differential method” was designed, in which one directly measures 
the wave functions using a position sensitive detector. In this method, a 
neutron guide of moderate height can be used and the deforming effect was 
kept rather small compared to the previous “integrated method”. The first 
attempt of the differential method was reported in [8], where a plastic nuclear 
track detector (CR39) with an uranium coating [9] was used as the detector. 

These quantum systems, in which the UCN bounces on a smooth floor, 
are expressed by the time-independent Schrodinger equation, 



( 1 ) 


with the linear potential 


V(z) 


mgz, z > 0, 
oo, z < 0, 


( 2 ) 
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where ijj n and E n are the eigenfunctions and eigenenergies, U is the reduced 
Planck constant, m is the neutron mass, and g is the standard gravitational 
acceleration. Equation (JT]) can be rewritten in dimensionless form as 

(i“£n)^n(£n)=°, (3) 

where = z/z 0 — E n /E 0 and the characteristic length and energy are given 
by 

2 » = ( 2^ )1/3 ~ 6 " m - (4) 

E 0 = (=£*)■/» ~0.6peV. (5) 

This is the Airy equation and solutions are described by Airy special func¬ 

tions, Ai(£ n ) and Bi(£ n ). Calculated probability distributions for the first five 
states are illustrated in Fig. |T| Note that = 0, where z n = z 0 E n /E 0 , are 
the classical turning points of a classical bouncing ball, which are denoted in 
the cross points of each eigenenergy and the potential line of mgz in Fig. [0 



Figure 1: Probability distributions of gravitationally bound states of neu¬ 
trons for the first five states Pi- Horizontal lines indicate eigenenegies E n 
and vertical lines show the classical turning points z n for each state. 

This system is a suitable device with which to test the inverse square 
law of standard Newtonian gravity and to search for new gravity-like short- 
range interactions around these scales. The first limit to a non-Newtonian 
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force with a Yukawa-type interaction potential using this microscopic system 
was reported in HU- The limits for a CP-violating Yukawa-type potential 
were shown in mm- Several experimental schemes have been proposed 
to improve the sensitivity to such hypothetical new physics. One idea to 
achieve a better resolution on the measurement of the characteristic length is 
to utilize a convex reflection mirror to magnify the neutron distributions m 
[T5] . In this review, we discuss an experiment using this scheme to precisely 
measure the UCN position-space wave functions, performed by the group of 
S. Komamiya of the University of Tokyo |16j . Experiments exploiting other 
ideas for measuring the energy scale, i.e. the energy differences between 
quantum states, are proposed by observing resonance transitions induced by 
a magnetic held [171120] and mechanical vibrations [2Tj. These projects are 
called GRANIT and qBounce, respectively. The first measurement of the 
resonance transition from the ground to the third state was reported by the 
group of H. Abele in 122 f H. For details of those projects, please see review 
papers [23] and [24] in this special issue. 

Precision measurement and issues related to this anal¬ 
ysis 

The precision measurement with a convex magnification mirror [16] was per¬ 
formed using a UCN source provided at ILL(Institut Laue-Langevin) [25] . 
The velocity distribution was measured by a TOF method and is well char¬ 
acterized by a Gaussian distribution with mean of 9.4 m/s and standard de¬ 
viation of 2.8 m/s. Figure^a) shows a schematic drawing of the experiment. 
The entire system was mounted on an anti-vibration table and magnetically 
shielded by a Permalloy sheet. The main components can be separated into 
three parts: a vertically thin collimating guide, a magnification mirror (rod), 
and a pixelated position sensitive detector. The layout of these components 
is shown in Fig. [2(b). 

Collimating guide 

The guide settles gravitationally bound quantum states. To clearly distin¬ 
guish each quantum state, which have energy differences of order of A E ~ 
1 peV, a resolving time of At ~ h/AE ~ 1 msec is required. Considering 
that the horizontal velocities of UCNs which we use are less than 15 m/s, the 
collimating guide was designed to have a sufficient length (192 mm in our 

1 Recently, a new limit for the CP-Violating Yukawa-type potential using the resonance 
method was reported in [25] . It also shows a limit for the chameleon field [27H251 , a dark 
energy candidates. 
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Figure 2: Schematic drawing of the precision measurement experiment using 
a convex magnification mirror pTO]. (a) is an overview of the system and (b) 
shows the geometry of the main components. 
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setup) to form each quantum states. The bottom mirror is made of polished 
glass with a roughness of R a (arithmetic average) = 0.03 /irn. The ceiling 
scatterer is a Gd-Ti-Zr alloy (54/35/11) deposited on glass, with a Fermi po¬ 
tential tuned to be nearly zero, and a roughness of R a = 0.4 /irn [30], which 
scatters out neutrons in the higher states. The guide selects lower states 
with appropriate populations to improve the contrast of the quantum spa¬ 
tial modulation^. The neutron loss models for the scatterer or rough surface 
have been discussed in detail in [3TM33] . and they are still interesting issues 
not only for UCN guiding applications but also for storage experiments such 
as neutron EDM measurements. In the experiment reviewed in this paper, 
empirical models of the loss rates are adopted [16]. The loss rate by the scat¬ 
terer, T n , is assumed to be proportional to the probability of finding neutron 
in the roughness region, and is given by 

r„ = 7 / _\ip n \ 2 * dz , (6) 

J h—25 


where 7 is a scaling constant, estimated from data to be 9.5ig;g x 10 4 s -1 , 
h is the height of the guide (100 /irn) , 5 is the roughness of the scatterer 
(0.4 /irn) , and ifj n are deformed wave functions in the guide. Neutron losses 
at the bottom mirror due to absorption, non-specular reflection, up scattering 
and other processes are modeled empirically as 


B 


n 




(7) 


in which the loss rate is assumed to be proportional to the classical bouncing 
number per unit time, g/2v Z)rLtmaxi where j3 is a scaling constant (estimated 

to be 0.38lo;o 4 ), v ZjrLtmax = \j2E n /m is the maximum vertical velocity of a 
neutron in the nth state, and E n is the eigenenergy of the deformed wave 
functions. The transmissivity of the guide for each state can be written as 


Pn OC 


exp[——(T n + B n )] 



( 8 ) 


where < ... > Vx indicate the average over the neutron horizontal velocities. 
By applying a diabatic transition from region II to HI, the population distri¬ 
bution (a probability in the paper [16]) of neutrons for each state at the end 
of region HI is estimated as in Fig. Ha). 

2 The eigenfunctions modulate somewhat coherently, especially around the lower region, 
except for the ground state, as seen in Fig. [T] To improve the contrast of the quantum 
modulation, a method using a negative step of several tens of /im, which transfers neutrons 

from the ground state to higher states, was considered in [5]. Results were reported 

in (TTIfTHI . 


6 






Principal quantum number n Height z(gm) 


Figure 3: (a) is an estimated population distribution (probability distribu¬ 
tion of the guide transmission) at the end of region IE. and (b) shows the 
corresponding Wigner phase space distribution [16]. 


Magnification mirror 

A Ni coated cylindrical rod is used as a magnifying convex mirror. Its radius 
is 3 mm and it is placed to have a grazing angle of 20 deg. for horizontally 
moving neutrons at the bottom floor level (z = 0). Figure [4] shows the 
magnification power as a function of the height [10] . ft gives about 20 times 
magnification around the lower region of z ~ 20 /irn. Before depositing the Ni 
coating, the cylinder was polished at the Research Center for Ultra-Precision 
Science and Technology, Osaka University. The roughness of the rod after 
the depositing was measured to be R a = 1.9 nm, two orders of magnitude 
smaller than the neutron wavelength. 

To calculate the detected position distribution on the detector surface, the 
wave function at the end of the bottom floor is re-expressed by the Wigner 
phase space distribution mm, 

roo £ 

W(z,p z ) = 2nhJ-oc 6ri e ~* PZV < " + 2 r, ^ Z ~ 2 V > ’ 

where p z is the momentum and p is a density operator. The Wigner distri¬ 
bution is known as a phase space formulation of quantum mechanics, and is 
widely used for quantum optics, for example in the study of decoherence [39] . 
As an application to massive particles, one can find a paper which shows a 
phase-space tomography of the Wigner distribution for a coherent atomic 
beam in a double-slit experiment [10]. Figure EI[b) shows the Wigner distri¬ 
bution constructed from the estimated populations in our experiment [16] . 
The time evolution of the Wigner distribution is calculated by the evolution 
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Figure 4: Magnification power M of the cylindrical rod as a function of the 
height ^ ra- It is about 20 times magnification around the lower region. 


of the density operator described by the Liouville - 

von Neumann equation 

dp i r tj - 

dt ’ 

(10) 

where 

~ u 2 

H = t + V 

By evaluating the kinetic part as 

(11) 

T = -=^-W(z,pO , 
m oz 

(12) 


and the potential part as 


oo oo 


£<* = £ 

1=0 1=0 


(-l) l (h/2) 21 d 2l+1 V(z) d 2l+l 
(2 1 + 1)! dz 2l+1 dpf +1 Z, ^ z ’ 


(13) 


one can obtain the Quantum Liouvillc Equation for the Wigner distribution 

m 




(14) 


In the case of V(z) = mgz, the right hand side of the equation vanishes 
and it becomes the classical Liouville equation. Therefore, in region IV, one 
can treat the evolution of each phase point of the Wigner distribution as a 
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classical path under gravity. Figure 0 shows the measured data and the best 
theoretical estimation using this model. The corresponding p-value is 0.715, 
and the experimental data support the quantum features described by the 
phase space formulation using the Wigner distribution [T6j . 



Figure 5: Black crosses show measured neutron position distribution [H)j . 
Modulated distribution was clearly measured. The result shows good agree¬ 
ment with the quantum expectation calculated using phase space formaliza¬ 
tion by Wigner distribution (p-value is 0.715). The best estimated line is 
shown by a red solid curve and the other lines indicate contributions of each 
state to the distribution. For definitions of fitting parameters, best fit values, 
systematic uncertainties, and the other details, see |16j. 


Position sensitive detector for UCNs 

A back-thinned CCD (HAMAMATSU S7030-1008) with thin Ti- 10 B-Ti layers 
is used for the position sensitive detector in the experiment [IBj. Its pixel 
size is 24 /irn x 24 /irn and the thickness of the active volume is about 20 /irn. 
The 10 B layer converts neutrons into charged particles by the nuclear reaction 
10 B(n,o;)'Li. The secondary particles are emitted in a nearly back-to-back 
configuration. One of them deposits its kinetic energy in the active area 
and creates a charge cluster, which typically spreads into nine pixels. The 
weighted center of the charge cluster is a good estimation of the incident 
neutron position. The thickness of the layers are 20 nm and 200 nm for Ti 
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and 10 B, respectively, and they are formed by evaporating directly on the 
CCD surface. The spatial resolution is measured to be 3.35 ± 0.09 //rri by 
evaluating the line spread function (LSF) (see Fig. [6]) using very cold neutron 
beams at ILL [TO] , For the detail of the evaluation scheme, see mm- A 
neutron converter of 6 Li using 6 Li(n,a) 3 H reaction was also investigated. It 
is concluded that the use of 10 B gives better spatial resolution £42]. 


Gd 


min 
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0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

Position (mm) 


Figure 6: Fitting results of a Gd shadow pattern [10]. Spatial resolution is 
measured to be 3.35 ±0.09 /ini by evaluating these line spread functions. For 
details of evaluation method, see [ 151142] . 

Other detectors using a silicon pixel device, Timepix [43], with 6 LiF and 
10 B converter were studied in [44]. Its pixel pitch is 55 /mi x 55 fin i and the 
thickness of its silicon layer is about 300 /mi. The spatial resolution for the 
6 LiF converter was evaluated by LSF to be 2.3 /mi, corresponding to about 
5.3 /mi in FWHM of the point spread function (PSF). The performance of 
the 10 B converter is estimated by Monte-Carlo simulations to be better than 
3 /mi in FWHM of the PSF. 

Another detector concept of Uranium coated plastic nuclear tracker (CR39) 
[9j was used in several experiments [8l l34ll36ll45] . Usually two fission frag¬ 
ments are emitted from a thin 235 U coating, and one of the daughter nuclei 
makes a track of defects in the CR39. By chemical etching, the diameter 
of track points is increased up to 1 /mi, allowing us to scan these vertexes 
using an optical microscope. Position resolution is around 1 /mi 0. By 
carefully analyzing the vertex shape, the spatial resolution can be improved 
to 0.7 /mi [45]. 
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Summary 

The quantum system of a gravitationally bound neutron is one of the most 
suitable tools to investigate gravity or gravity-like hypothetical interactions 
around the scale of 10 /an in length or 1 peV in energy. After establishing 
this research field by the pioneering work in observation of the quantum state 

[6], experimental schemes for precision measurements of these characteristic 
scales have developed rapidly, and nowadays, one can establish limits for 
parameter spaces of new physics experimentally [TTHT3[|26]. Furthermore one 
can start discussing a phase space formulation of quantum physics for the 
gravitationally bound quantum state [T6] . An experiment for a possible test 
of phase space formalization using spatial interference is under preparation. 

Acknowledgments 

The authors would like to thank Valery V. Nesvizhevsky (Institut Laue- 
Langevin), Hartmut Abele (Vienna University of Technology), William Snow 
(Indiana University), Peter Geltenbolt (Institut Laue-Langevin), and all par¬ 
ticipants in the GRANIT-2014 Workshop for interesting discussions and help¬ 
ful suggestions. This material is based upon work supported by JSPS KAK- 
ENHI Grants No. 20340050 and No. 24340045, and Grant-in-Aid for JSPS 
Fellows No. 22.1661. 


References 

[1] R. Colella, A. W. Overhauser, and S. A Werner, Phys. Rev. Lett. 34, 
1472 (1975) 

[2] J. J. Sakurai, Phys. Rev. D 21, 2993 (1980) 

[3] J. L. Staudenmann et al ., Phys. Rev. A 21, 1419 (1980) 

[4] S. A. Werner et al, Physica B 151 , 22 (1988) 

[5] K. C. Littrell, B. E. Allman, and S. A. Werner, Phys. Rev. A 56 , 1767 
(1997) 

[6] V. V. Nesvizhevsky et al, Nature 415 , 297 (2002) 

[7] V. V. Nesvizhevsky et al, Phys. Rev. D 67 , 102002 (2003) 

[8] V. V. Nesvizhevsky et al, Eur. Phys. J. C 40 , 479 (2005) 


11 


[9] V. V. Nesvizhevsky et al, Nucl. Instrum. Meth. A 440, 754 (2000) 

[10] G. Ichikawa, PhD. Dissertation, The University of Tokyo (2013) 

[11] V. V. Nesvizhevsky and K. V. Protasov, Class. Quantum Grav. 21, 4557 
(2004) 

[12] S. BaeBler it et ah, Phys. Rev. D 75, 075006 (2007) 

[13] S. BaeBler it et ah, Nucl. Instrum. Meth. A 611, 149 (2009) 

[14] V. V. Nesvizhevsky, Preprint ILL 96NE14T (1994) 

[15] T. Sanuki et al. , Nucl. Instrum. Meth. A 600, 657 (2009) 

[16] G. Ichikawa et al. , Phys. Rev. Lett. 112, 071101 (2014) 

[17] P. Schmidt-Wellenburg et al. , Nucl. Instrum. Meth. A 611, 267 (2009) 

[18] M. Kreuz et al, Nucl. Instrum. Meth. A 611, 326 (2009) 

[19] S. BaeBler et al. (GRANIT Collaboration), C. R. Physique 12, 707 (2011) 

[20] V. V. Nesvizhevsky (GRANIT Collaboration), Mod. Phys. Lett. A 27, 
1230006 (2012) 

[21] H. Abele et al. , Phys. Rev. D 81, 065019 (2010) 

[22] T. Jenke et al. , Nature Phys. 7, 468 (2011) 

[23] D. Roulier et al. , Adv. High En. Phys. this issue (2014) 

[24] G. Cronenberg et al. , Adv. High En. Phys. this issue (2014) 

[25] A. Steyeri et al, Phys. Lett. A 116, 347 (1986) 

[26] T. Jenke et al., Phys. Rev. Lett. 112, 151105 (2014) 

[27] J. Khoury and A. Weltman, Phys. Rev. Lett. 93, 171104 (2004) 

[28] D. F. Mota and D. J. A. Shaw, Phys. Rev. Lett. 97, 151102 (2006) 

[29] J. Khoury, Class. Quantum Grav. 30, 214004 (2013) 

[30] Y. Kamiya et al, KURRI Progress Report COl-1 (2008) 

[31] A. E. Meyerovich and V. V. Nesvizhevsky, Phys. Rev. A 73, 063616 
(2006) 


12 


[32] A. Yu. Voronin et al., Phys. Rev. D 73, 044029 (2006) 

[33] A. Westphal et al., Eur. Phys. J. C 51, 367 (2007) 

[34] V. V. Nesvizhevsky et al, ILL Annual Report (2004) 

[35] H. Abele et al, Nucl. Phys. A 827, 593c (2009) 

[36] T. Jenke et al, Nucl. Instrum. Meth. A 611, 318 (2009) 

[37] E. Wigner, Phys. Rev 40, 749 (1932) 

[38] M. Hillery et al., Phys. Rep. 106, 121 (1984) 

[39] L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cam¬ 
bridge University Press 1995) 

[40] Ch. Kurtsiefer, T Pfau, and J. Mlynek, Nature 386, 150 (1997) 

[41] W. P. Schleich, Quantum Optics in Phase Space (Wiley-VCH 2001) 

[42] S. Kawasaki et al, Nucl. Instrum. Meth. A 615, 42 (2010) 

[43] X. Llopart et al, Nucl. Instrum. Meth. A 581, 485 (2007) 

[44] J. Jakubek et al, Nucl. Instrum. Meth. A 600, 651 (2009) 

[45] T. Jenke et al, Nucl. Instrum. Meth. A 732, 1 (2013) 


13 


